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We propose and analize a graphene tunneling transit time device based on a heterostructure with 
a lateral p-i-n junction electrically induced in the graphene layer by the applied gate voltages ol 
different polarity. The depleted i-section ol the graphene layer (between the gates) serves as both 
the tunneling injector and the transit region. Using the developed device model, we demonstrate 
that the ballistic transit ol electrons and holes generated due to interband tunneling in the i-section 
results in the negative ac conductance in the terahertz frequency range, so that the device can serve 
as a terahertz oscillator. 



Graphene is considered as a promising candidate for 
different future electronic and optoelectronic devices. Its 
most distinctive features beneficial for device applica- 
tions include high electron and hole mobilities in a wide 
range of temperatures, the possibility of bandgap engi- 
neering (creation of the graphene-based structures with 
the energy gap from zero to fairly large values), forma- 
tion of electrically induced lateral p-n junctions (see, for 
instance, 0, 0, 1 S i ) . 

The operation of graphene field-effect transistors (G- 
FETs) is accompanied by the formation of the lateral n- 
p-n (or p-n-p) junction under the controlling (top) gate 
and the pertinent energy barrier 

HUSH- The tunnel- 
ing across such a n-p-n junction prevents the achievement 
of a low off-state current @ . This limits possible realiza- 
tion of G-FETs in large scale digital electronic circuits 
and forces to consider the graphene structures in which 
the energy gap is reinstated (graphene nanoribbons and 
graphene bilayers with the ene rgy gap open by the trans- 
verse electric field) [13, El [H UMLM, EH Ell • 

Unique properties of graphene already produced not 
only using peeling technology but also epitaxial methods 
as well as graphene nanoribbons and bilayers, particu- 
larly, experimental evidences of the possibility of ballis- 
tic electron and hole transport in samples with several 
micrometer sizes even at room temperatures (see, for in- 
stance, Refs. pj], EII) stimulate inventions of different 
graphene-based devices which could not be realized the 
past using the customary materials. 

In this paper, we propose a transit-time oscillator 
which can operate in the terahertz (THz) range of fre- 
quencies and substantiate the operational principle of the 
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device. The operation of the device in question is asso- 
ciated with the tunneling electron injection in an elec- 
trically induced reverse biased lateral p-i-n junction and 
the electron and hole transit-time effects in its depleted 
section. In the following, this device is referred to as 
the graphene tunneling transit-time (G-TUNNETT) ter- 
ahertz oscillator. The tunneling generation through the 
zero energy gap and propagation of electrons and holes 
with their directed velocity v x close to the characteristic 
velocity vyy ~ 10 8 cm/s of the graphene energy spec- 
trum (v x ~ vw), can provides significant advantages of 
G-TUNNETTs in comparison with the existing and dis- 
cussed TUNNETTs based on the conventional semicon- 
ductor materials (see, for instance, Refs. [H, [2(| and the 
references therein). Using the developed device model, 
we calculate the high-frequency characteristics. The de- 
vice under consideration comprises a graphene layer with 
the source and drain contacts and the gates. The lateral 
p-i-n junction can be formed by two top gates biased by 
the dc voltages V p and V n of different polarities (V p < 0, 
V n > 0). The G-TUNNETT structure under considera- 
tion and its band diagram (at the applied source drain- 
voltage V, so that the lateral p-i-n junction is reverse 
biased) are schematically shown in Fig. la and Fig. lb, 
respectively. It is assumed that apart from a dc compo- 
nent Vq > which corresponds to a recerse bias, the net 
source-drain voltage V comprises also an ac component 
SV exp(— iut): V = Vq + SV exp(— iu t), where SV and 
w are the signal amplitude and frequency, respectively. 

Thus, the graphene layer is partitioned into three sec- 
tions: p- and n- sections adjacent to the source and drain 
contacts, respectively, and depleted section in the center 
(i-section). The depleted section of the graphene layer 
plays the dual role: (1) a strong electric field in the 
i-section provides the tunneling generation of electrons 
and holes and (2) this section serves as the transit re- 
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FIG. 1: Schematic view of a G-TUNNETT structure (a) and 
its band diagram (b) . 



gion where the generated electrons and holes propagating 
ballistically across this region induce the current in the 
p- and n- sections (as well as, possibly in the contacts) 
and, hence, the terminal current in the source-drain cir- 
cuit. Considering for definiteness the "symmetric" device 
structure with V n — — V p — V t ) corresponding to Fig. 1, 
we shall calculate its source-drain ac conductance a^f 
(dynamic conductance) as a function of the structural pa- 
rameters and the signal frequency and demonstrate that 
crff < in certain (THz) frequency ranges. 

Usingthe value of the interband tunneling probabil- 
ity [f| [2l[, and roughly estimating the electric field in 
the i-section as £ ~ V/21, where 21 is the length of the 
depleted i-section (21 ~ D if the spacing between the top 
gates D 3> Wti were Wt is the thickness of the layers 
separating the graphene layer and the gates; more strict 
and detailed calculations can be found in Refs. [2~il. |22|). 
for the rate of the tunneling generation of electron hole 
pair (per unit length in the transverse direction) in the 
i-section in the situation under consideration one can ob- 
tain 
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Here go = 2e 2 /irh, e is the electron charge, and h is 
the reduced Planck constant. This corresponds to the 
source-drain dc current 



J = 2eG = 
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The ac component of the generation rate is given by 
< n 3ffo ta 1 / 2 -^ 3 J ° 
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The tunneling probability w is a fairly sharp function 
of the angle between the direction of the electron (hole) 
motion and the x-direction (from the source to the drain) : 
w{9) ~ exp(— asm 2 9), where a oc D is rather large [f|. 
Considering this, one can disregard some spread in the 
x-component of the velocity (A v x /v\y — l/2a) of the in- 
jected electrons and assume that all the generated elec- 
trons and holes propagate in the x-direction with the 
velocity v x ~ vw- 

In the case of ballistic electron and hole transport in 
the i-section when the generated electrons and holes do 
not change the directions of their propagation, the conti- 
nuity equations governing the ac components of the elec- 
trons and holes densities 5Ej and can be presented 
as: 
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The boundary conditions are as follows: 8T,^\ x=T i = 0. 
Solving Eq. (4) with these boundary conditions, we arrive 
at 



8G^ exp[iuj(l ± x)/vw] — 1 
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Considering this, for the net ac current 5 J u (x) = 
ev w [S^Z( X ) + ^i( x )]> created by both the generated 
electrons and holes, one can obtain 



5J u (x)=e6Gl 



e tuJT cos{uTx/l) - 1 
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Here, we have introduced the characteristic transit time 
t = l/vw Since the real time of the electron and hole 
transit across the i-section varies from to 21 /vw, t is 
actually the mean transit time. If the signal frequency 
lo is smaller than the plasma oscillations in the p- and 
n-sections f2, (uj < ft), these sections can be considered 
as just highly conducting electrodes, so that the terminal 
ac current is mainly induced in them. In this case, the ac 
component of the source-drain terminal current SJ^f, i.e., 
the ac current induced by the propagating electrons in 
the external circuit connecting the source and the drain 
can, according to the Ramo-Shockley theorem [23|,[24[, be 

presented as SJ^f — f_ t dx g(x)5J u} (x) — iu; C8V Ul where 
g(x) is the form-factor determined by the shape of highly 
conducting regions (the p-section and the drain contact) 
and C is the i-section geometrical capacitance. The fac- 
tor g(x) is associated with different contributions to the 
induced current of the electrons and holes at different dis- 
tances from the highly conducting regions. For the bulky 
conducting regions g(x) ~ 1/2/, whereas for the blade- 
like conducting regions [III, g(x) 
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the i-section capacitance one can use the following for- 
mula: [26] C = (as/47r 2 ) ln(4L/Z), where as is the dielectric 
constant. Considering Eqs. (3), we arrive at the following 
formula for the ac conductance a^f = 8J^ d /SV u : 
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FIG. 2: Real part and imaginary part (inset) of the ac con- 
ductance alf as functions of signal frequency / = cj/2tv for 
different lengths of i-section (different transit times). 
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(x = 0.25 ps, c = 0.34) 
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(x = 0.35 ps, c = 0.30) 
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where <tq = Jq/Vq is the dc conductance and J7b(£) is the 
Bcssel function. 

The real and imaginary parts of the ac conductance 
are given by 
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where c = 2C/3ero T - At relatively low signal frequencies 
(wt < 1), Eqs. (7) and (8) yield Re erf, d ~ (3cr /2)(l - 
7w 2 r 2 /12) and Im a s J ~ (3<7n)(3/4- c)wr, respectively. 
Thus, (Jo d — 3cro/2, i.e., is equal to the differential dc 
conductance cr^ftf = {dJo/dVo)\vt=const- Assuming that 
21 = 0.5 fim, 2L = 2 fim, ae = 3.8, and V = 0.1-1 V, 
we obtain er - (2.09 - 6.64) x 10 12 s" 1 and c ~ 0.34. 
The real part of the ac conductance Re aff turns zero at 
the following signal frequencies / = uj/2tt: ~ 0.38/r, 
ft = 0.5/r, ~ 0.88/t, /+ = 1/t,... In the frequency 
ranges 0.38/r < / < 0.5/r and 0.88/r < / < 1/r, the 
real part of the ac conductance is negative. In partic- 
ular, at 21 = 0.7 urn (r = 0.35 ps), /f ~ 1.08 THz 
and fy = 1.43 THz. The real part of the ac con- 
ductance Reerf, d reaches minima at certain frequencies 
/i,/2,..., which fall into the intervals /-f < /i < /{*", 
< fi < f^v The quantity Rea^f increases with in- 
creasing source-drain voltage and, consequently, dc cur- 
rent. A decrease in the i-section length at fixed source- 
drain voltage gives rise to an increase both in |Recrf, d | and 
in the frequencies /i,/2,.-., where Recr^ d exhibits man- 
ima. 



FIG. 3: Real part vs imaginary part dependence with lot as 
a parameter for different length of i-section (different transit 
times and normalized capacitances). 



Figure 2 shows the real and imaginary parts of the 
ac conductance a^f as functions of the signal frequency 
/ = lj/27t calculated for 21 — 0.5 /im and 21 = 0.7 /an. 
The real part exhibits a pronounced oscillatory behavior 
with the frequency ranges where it has different signs. 
The frequency dependence of the imaginary part corre- 
sponds to the domination of the inductive component 
at low frequencies (associated with the contribution of 
electrons and holes) and the capacitive component at el- 
evated frequencies (due to the contribution of the geo- 
metrical capacitance) . The relative value of the real part 
of the ac conductance at the first minima at 21 = 0.5 /urn 
is fairly moderate: R-eo^L^ * /ctq ~ —0.034. However, 
it is markedly larger than that in some new concept 
THz devices [27]] • Considering the above estimate of 
00 at Vo = 0.1 — 1 V, one can obtain Re£7 , * < L 2 7r A — 
-(0.71 - 2.26) x 10 11 s _1 . As follows from the above 
estimates and Fig. 2, the negativity of the real part of 
the ac conductance at the frequencies / > 1 THz can 
be achieved in the G-TUNNETT structures with the i- 
scction length only moderately smaller than one microm- 
eter. We assumed above that the generated electrons 
and holes in the i-section propagate ballistically. Taking 
into account that the electron-electron and electron-hole 
scattering in this (depleted) section should be weak, the 
mean-free paths in graphene can be of several microme- 
ters even at room temperatures [l7], [U[ . Our estimates 
show that the effect of the electron and hole space charge 
on the potential distribution in the i-section is weak and 
the pertinent plasma phenomena arc insignificant. 

The fundamental plasma frequency in the gated p- and 
n-sections, in which the spectrum of the plasma waves is 
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a sound- like f2 = ws/2{L — I), where s is the charac- 
teristic velocity, which at sufficiently high gate voltages 
when the electron (hole) Fermi energy in the n-section (p- 
section) Ef — (fi v w/2V2) y/ssVt /eW t ksT, is given 
by s ~ y / 2e 2 W t eF/eeh 2 > vw [28], 2L is the spacing 
between the source and drain contacts, T is the temper- 
ature, and ks is the Boltzmann constant. Since the the 
plasma frequency in graphene can be rather large (owing 
to a large s) , the condition u> < fl assumed above can be 
easily fulfilled for the signal frequencies in the THz range 
(as well as the condition u> ~ O) by a proper choice of 
the device structure length 2L and the gate voltage Vt- 
Indeed, setting, for example, s = (3 — 6) x 10 8 cm/s, 
21 = 0.5 //m, and 2L = 1.5 — 2 fim, one can obtain 
fl/2TT =1-3 THz. 

If /f,/+ < Q/2n (or /„-,/+ < 0/2tt), a G- 
TUNNETT coupled with the proper resonant cavity can 
serve as a THz oscillator. However, when /£" < n/2ir < 
fi (or / 2 ~ < ri/27r < ); i- c -; when the frequency of the 
plasma oscillations falls into the range where Re a^f < 0, 
the hole and electron systems in the gated p- and n- 
sections, respectively, can play the role of the resonators. 
In this situation, the self-excitation of the plasma oscil- 
lations (plasma instability) can be possible if the quality 
factor of these oscillations Q ~ Fl/v, where v is the fre- 
quency of hole and electron collisions with defects and 
acous tic p honons in the p- and n-sections, is sufficiently 
large [26| (see also Refs. [HI [HJ). It should be noted 
that the electron-electron and hole-hole scattering pro- 



cesses do not affect Q except as via relatively weak ef- 
fects of electron and hole viscosity [3l|. As a result, in 
sufficiently perfect graphene layers, v can be small not 
only in the i-section but in the p- and n-sections as well. 
Thus, a G-TUNNETT can work as a THz source in the 
regime of plasma instability as well. In such a case, a G- 
TUNNETT should be supplied by an antenna, although 
the gates can also play the role of the latter. 

In conclusion, we have proposed a G-TUNNETT and 
calculated its ac conductance as a function of the sig- 
nal frequency and the structural parameter using the de- 
veloped device model. We have demonstrated that the 
ac conductance exhibits the frequency ranges where it is 
negative. Due to high directed energy independent ve- 
locities of the electrons and holes generated owing to the 
interband tunneling, these frequency regions correspond 
to the THz range at relatively large length of the i-region. 
A G-TUNNETT can work as an active element of THz 
oscillators with a complemetary resonant cavity or can 
immediately emit THz radiation (in the plasma instabil- 
ity mode). 
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